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Abstract 

Lemaitre - Tolman (L-T) toy models with a central observer have 
been used to study the effect of large scale inhomogeneities on the SN la 
dimming. Claims that a giant void is mandatory to explain away dark en- 
ergy in this framework are currently dominating. Our aim is to show that 
L-T models exist that reproduce a few features of the ACDM model, but 
do not contain the giant cosmic void. We propose to use two sets of data 
- the angular diameter distance together with the redshift-space mass- 
density and the angular diameter distance together with the expansion 
rate - both defined on the past null cone as functions of the redshift. We 
assume that these functions are of the same form as in the ACDM model. 
Using the Mustapha-Hellaby-Ellis algorithm, we numerically transform 
these initial data into the usual two L-T arbitrary functions and solve the 
evolution equation to calculate the mass distribution in spacetime. For 
both models, we find that the current density profile does not exhibit a 
giant void, but rather a giant hump. However, this hump is not directly 
observable, since it is in a spacelike relation to a present observer. The 
alleged existence of the giant void was a consequence of the L-T mod- 
els used earlier because their generality was limited a priori by needless 
simplifying assumptions, like, for example, the bang-time function being 
constant. Instead, one can feed any mass distribution or expansion rate 
history on the past light cone as initial data to the L-T evolution equation. 
When a fully general L-T metric is used, the giant void is not implied. 
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1 The historical background of the problem 



In the framework of the homogeneous and isotropic standard cosmological model, 
the dimming of the type la supernovae as compared to their expected luminosity 
in an Einstein dc Sitter model is interpreted as a consequence of an assumed 
accelerated expansion of the Universe. This leads to the widespread belief in a 
'dark energy' component currently dominating the energy budget of our Uni- 
verse. But this is not the only possible explanation of the SN la observations. 

Shortly after the discovery by Riess et al. (1998) and Perlmutter et al. 
(1999) of the supernova dimming, it was proposed by several authors that this 
effect could be due to the large-scale inhomogeneities (Pascual-Sanchez 1999; 
Celerier 2000; Tomita 2000; 2001a & b). After a period of relative disaffection, 
this proposal experienced a renewed interest about five years ago. 

Three methods have been used to implement such a proposal: computation 
of backreaction terms in the dynamical equations using an averaging procedure 
proposed by Buchert (2000, 2001), calculations in the framework of a perturba- 
tive scheme and the use of exact inhomogeneous models, in particular those of 
Lemaitre (1933) - Tolman (1934) (L-T) (see Celerier 2007, for a review). 

The L-T model became rapidly popular for the purpose of mimicking 'dark 
energy' because it exhibits three interesting features: i) it is one of the few exact 
solutions of General Relativity able to represent a physically consistent model 
of the matter dominated era of the Universe, ii) among these few, it is the most 
easily tractable from a computational point of view, iii) it is not an alternative 
to, but a generalisation of the Friedmann dust models (which are contained in it 
as a subcase), so can reproduce all the Friedmann-based results, including those 
of the "concordance" ACDM model, with an arbitrary precision. For more on 
this, in relation to the main subject of this paper, see the last section. 

Three classes of models have been constructed with the L-T solution: i) 
models where the observer is located at the centre of a single L-T universe 
(e.g., Iguchi et al. 2002; Alnes et al. 2006; Apostolopoulos et al. 2006; Bolejko 
2008; Garcia-Bellido & Haugb0lle 2008a, 2009), ii) models where the observer is 
located off the centre of such a universe (e.g., Schneider & Celerier 1999; Apos- 
tolopoulos et al. 2006; Alnes & Amarzguioui 2007) hi) Swiss-cheese models 
where the holes are L-T bubbles carved out of a Friedmannian homogeneous 
background (e.g., Brouzakis et al. 2007, 2008; Biswas & Notari 2008; Marra et 
al. 2007). 

As will be recalled in Sec. [H an L-T model is defined by two independent 
arbitrary functions of the radial coordinate, which can be fitted to the observa- 
tional data. However, in most of the models currently available in the literature, 
the authors have artificially limited the generality by giving the L-T initial-data 
functions a handpicked algebraic form (depending on the authors' feelings about 
which kind of model would best represent our Universe), with only a few con- 
stant parameters being left arbitrary - to be adapted to the observations. 

Another way in which the generality of the L-T models was artificially lim- 
ited was the assumption that the age of the Universe is everywhere the same, 
i.e. that the L-T bang-time function ts is constant. With ts being constant, 
the only single-patch L-T model that fits observations is one with a giant void 
(Iguchi et al. 2002; Yoo et al. 2008). Conceptually there is nothing wrong with 
a non-simultaneous big bang (even though this is a radical qualitative differ- 
ence with the FLRW models), but one should exercise caution when referring 
to (t — £b) as the actual age of the Universe. The L-T model is too simple to 

1 However, these authors have shown that the CMB data put very stringent limits on 
the distance of the observer from the centre of the model or/and on the amplitude of the 
homogeneities in an off-centre observer model. 
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extend it up to instants earlier than decoupling. Therefore ts should merely be 
regarded as a function that describes a degree of inhomogeneity of the initial 
conditions rather than as the actual instant of birth of the Universe. 

The argument brought in defense of the constant ts assumption is this: a 
non-constant ts generates decreasing modes of perturbation of the metric (Silk, 
1977, Plebahski and Krasihski, 2006), so any substantial inhomogeneity at the 
present time stemming from ts, r 7^ would imply 'huge' perturbations of ho- 
mogeneity at the last scattering. This, in turn, would contradict the CMB 
observations and the implications of inflationary models (we deliberately do not 
give references here, to avoid blaming any single individual for what seems to 
be a piece of conventional wisdom). However, these are only expectations that 
should not be treated as objective truth until they are verified by calculations. 
Such calculations have already been done, and it turned out that the inhomo- 
geneities in tg needed to explain the formation of galaxy clusters and voids are 
of the order of a few hundred years (Bolcjko et al. 2005; Krasinski and Hellaby 
2004; Bolejko 2009). Then, on the basis of Bolejko's (2009) models 4 and 5, one 
can calculate that for a structure of present radius 30 Mpc this age difference 
between the oldest and youngest region would generate CMB temperature fluc- 
tuations equal to AT/T = 3.44 x 10~ 6 and AT/T = -2.35 x 10~ 6 , respectively. 
This is well-hidden in the observational errors at the present level of precision. 
(In the future, when, presumably, the precision will improve, these results may 
possibly be used to measure the gradient of ts) So there is no observational 
justification to the assumption ts = constant. 

These are the reasons why, in the recent years, we have seen the increase in 
popularity of void models, where the observer is located at or near the centre 
of a large, huge, giant L-T void of size of up to a few Gpc. Many authors have 
constructed classes of L-T models with a central local void and have shown that 
they were able to fit the SN la and other cosmological data provided the void is 
large (e.g., it has a diameter of 400 h^ 1 Mpc in Alexander et al. 2007) or even 
huge (e.g., 1.35 Gpc in Alnes et al. 2006, > 2 Gpc in Garcia-Bellido & Haugb0lle 
2008a and 2009), depending on the features of the particular model they had 
chosen. This contributed to the spreading of the belief in the necessity of a 'giant 
local void' to resolve the 'cosmological constant problem' with L-T models. 
However, as shown by Mustapha et al. (1998) and used as an illustration for 
the application to the supernova data and the 'cosmological constant problem' 
by Celericr (2000), a given set of isotropic data can constrain only one of the 
two free functions of an L-T model and therefore, after fitting the supernova 
data with a given L-T solution, we are left with plenty of room to accommodate 
more observations. 

Actually, a few authors discarded the central void hypothesis and proposed 
models with no such void (e. g., Iguchi et al. 2002). Enqvist and Mattsson 
(2007) even showed that the fitting of the SN la data can be better with L-T 
models where the density distribution is constant on a constant-time hypersur- 
face than with the ACDM model (see also Bolejko 2008 and Bolcjko & Wyithe 
2009). Even though such a density distribution is not what is actually observed 
at very large scales by astronomers (it is not even observable, being in a space- 
like relation to the central observer), the p(to,r) = const configuration vividly 
illustrates how misleading the FLRW-based geometrical intuitions can be. A 
spatial distribution of matter can radically change with time in consequence 
of an inhomogeneous expansion distribution in space. Our models will provide 
more examples of this phenomenon, and we will come back to this point in the 
conclusion. 

Our aim here is to show, using two explicit examples reproducing the obser- 
vational features of the ACDM model, that a giant void is not at all a necessary 
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implication of using L-T models. We propose to use input functions that can 
be derived from observations @. Our L-T toy models will be constrained by the 
angular diameter distance together with the redshift-space mass-density or the 
angular diameter distance together with the expansion rate. 

It should be noted that these functions have not the same form in the ACDM 
model and in giant-void L-T models. For example, let us consider the giant 
void model from Bolejko & Wyithe (2009) with radius of 2.96 Gpc and density 
contrast of 4.05. The redshift-space mass-density for this model and for the 
ACDM model are shown in Fig. [TJ As seen, at z ~ 1 the difference between 
these two models is more than a factor of 2. Also, the expansion rate as a 
function of redshift behaves differently (for details and constraints coming from 
H(z) see Bolejko & Wyithe, 2009). Thus giant void models have difficulties to 
mimic all the observational features of the dark energy model (Zibin et al. 2008; 
Clifton et al. 2009). 

In this paper, we show that if the observational data are properly fitted to 
these ACDM functions, then a giant void is not mandatory to explain them. 
In fact the L-T models that mimic our choice of observational features of the 
ACDM model have a central Gpc-scale overdensity rather than an underdensity. 
We emphasise that what we reproduce in our L-T model are not the actual 
observational relations, but the parameters of the ACDM model fitted to the 
observations - which is not the same thing. 

Note that this model is not designed to reproduce all the available cosmo- 
logical data, nor is it to be considered as the final model of our Universe. Its 
purpose is to exemplify the proper use of L-T models and to show what can 
come out of it. Moreover, it should be understood as tentative beyond the 
redshift range in which the ACDM functions we use are robustly established. 

As is usual in the study of L-T models, we choose to use a comoving and 
synchronous coordinate system for the majority of this work. Such a coordinate 
system is uniquely defined by the flow lines of the fluid and allows the line- 
element to be written in a simple form. However, it is of course the case that 
quantities such as energy density profiles on space-like volumes are sensitive 
to the choice of hypersurface on which they are recorded. To illustrate this 
dependence, and the effect of considering other foliations, we also present our 
final results on a set of hypersurfaces in which each fluid element is the same 
distance from the initial singularity along the world-lines of the dust particles. 
Such a choice allows us to consider the energy density of different regions when 
they are at the same age, and becomes a comoving and synchronous coordinate 
system in the constant bang time models where giant voids have often been 
inferred. 

In Sec.[5]we recall briefly the main properties of the L-T solution and give the 
equations to be integrated on the light cone. In Sees. 13.11 and 13.21 we describe 
both models and give the results of our numerical calculations. Section [3] is 
devoted to a discussion and a summary of our conclusions. 

2 One should be aware that there is a great deal of phenomenology involved in interpreting 
the observations. For example, with supernovae, the cosmological model predicts D^, while 
what is actually observed is the flux. We can deduce the absolute luminosity only on the basis 
of some empirical methods. Similarly with galaxy number counts - the cosmological model 
predicts m(z)n(z) where m(z) is an average mass per source and n(z) is number counts. The 
whole information about the galaxy evolution and their mergers is encoded in m(z) - however 
in galaxy redshift surveys we observe only n(z). In this paper we do not focus on the problem 
of observations and assume Dr, and ran as in the ACDM. 
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Figure 1: The redshift-space mass-density in the ACDM model and in a giant 
void L-T model (with radius of 2.96 Gpc and density contrast of 4.05). At the 
redshift around 1 one should observe twice as many objects as in the ACDM 
model. 



2 The Lemaitre — Tolman solution 

The L-T model is one of two classes of spherically symmetric solution^ of 
Einstein's equations where the gravitational source is dust. In comoving and 
synchronous coordinates, its line element reads 

ds 2 =c 2 dt 2 ^4 T dr 2 -i? 2 (t,r)(d l 9 2 +sin 2 z9d^ 2 ), (1) 

l + 2E(r) y ' r ' J 

where E(r) is an arbitrary function and i?, r = dR/ dr. With a vanishing cosmo- 
logical constant, which is the case we consider here, R(t,r) obeys the following 
first integral of one of the Einstein equations: 

where R,t = dR/dt and M = M(r) is another arbitrary function of integration. 
M(r) is the gravitational mass contained within the comoving spherical shell 
at any given r, while E(r) is the energy per unit mass of the particles on that 
shell. Moreover, E determines the space curvature at each r-value. The mass 
density in energy units follows from the other nontrivial Einstein equation and 
is 

2 2M, r 8nG 
KP ° = R 2 R ' 16 K = ~c^' ^ ' 

The solutions of J2]| can be written as 

R(t,r) = ^j-t(t,r), C ( t - tfl ( r )) = ^^_ €(t , r ), (4) 

where 



3 For the presentation of the other class see Plebanski and Krasinski (2006). It is called 
there the Datt - Ruban solution (Ruban 1968, 1969). It has interesting geometrical and 
physical properties, but so far has found no astrophysical application. 
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1. when E < (elliptic evolution): 



X (r) = -2E(r), 
4> = 1 — cos?7, 
£ = r]-smr], 



(5) 



2. when E = (parabolic evolution): 



X(r) = 1, 

e = We, 



(6) 



3. and when E > (hyperbolic evolution): 




(7) 



where ij is a parameter (dependent on i and r). The is(r) is another arbitrary 
function that appears as an integration 'constant' and is interpreted as the bang 
time, i.e., the Big Bang is not simultaneous for all values of the r coordinate. 

Since all the formulae given so far are covariant under coordinate transfor- 
mations of the form r = g(r), one of the functions E(r), M(r) and ts(r) can 
be fixed at will by the choice of g. Therefore, once this choice is done, a given 
L-T model is fully determined by two of these arbitrary functions. 

However, as shown by Mustapha et al. (1998), and used as an illustration 
for application to the supernova data and the 'cosmological constant problem' 
by Celerier (2000), a set of isotropic data corresponding to a given observable 
can constrain only one of the two free functions, and therefore the fitting of 
the supernova data, i.e. of the function Dl(z), with a given L-T solution, still 
leaves the other function free - and available for fitting to another set of data. 
Thus, we must also assume another set of initial conditions, e.g. the redshift- 
space mass-density m(z)n{z) or the expansion rate H(z). We will take these 
functions to be identical to the corresponding functions in the ACDM model, 
assuming they reflect the observational data. By this, we want to show that 
there is no antagonism between the inhomogeneous cosmology and the ACDM 
model, and that the first can predict the same results as the second even if 
A = 0. However, whenever possible we present the real data to show that there 
is still much room within observational errors for different profiles. 

Using the reciprocity theorem (Etherington 1933, Ellis 1971) the luminosity 
distance can be converted to the angular diameter distance 



For a ray issued from a radiating source and proceeding towards the central 
observer on a radial null geodesic the following equation holds 



and the equation for the redshift reads (Bondi 1947; Plebahski & Krasihski 
2006): 



cdi(r) 
dr 



R,r (t(r),r) 

V 1 + *W) 



(9) 
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dz(r) _lR,tr(t(r),r) 



l + z{r) dr c y /\ + 2E(r)' 

For later reference we will need also the following equation, which follows easily 
from © and (ITU1) : 

1 dz(t) _ R, tr (t,r(t)) 
l + z(t) dt R, r (t,r(t))' [ ' 



3 The L— T model with no central void 

Our aim is now to show we can design an L-T model with no central void able 
to reproduce the angular diameter distance-redshift relation as inferred from 
the SN la data, smoothed out as in the framework of a ACDM model. For this 
purpose, we propose to use the following additional conditions to specify the 
arbitrary functions M(r), E(r) and ts(r). 



3.1 The model defined by Da{z) and m(z)n(z) 
3.1.1 The MHE procedure 

The algorithm used to find the L-T model reproducing the Da(z) and n(z) 
data was first developed by Mustapha, Hellaby & Ellis (1997). Let us recall its 
major steps and equations. 

The radial coordinate r is chosen so that, on the past light cone of (t, r) = 
(*o,0), _ 

- R,r = 1. (12) 
VT+2E 

(Note: this choice of r is possible only on a single light cone. In the following, 
we always refer to the light cone of (t,r) = (to,0).) This choice of coordinates 
simplifies the null geodesic equation: 



ct(r) = do — r, 



(13) 



where we denote quantities on this null cone by a hat. Furthermore, ([3]) now 
becomes 



{nc 2 /2)pR 2 = 



Vl + 2£ 

The total derivative of the areal radius R gives 

dR -g- g- dt 

— — — R, r + R,t 3~ ■ 

dr dr 
Using ©, (H2) and (fT5)) . the above equation can be written as 



(14) 



(15) 



dR 

dr 



VT+2E = — R, t =T\ 



I2M 
R 



2E. 



(16) 



This can be solved for E(r) 



1 + 2E=<- 



dR 
dr 




(17) 
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Using (fT4"|) the above becomes 

KC 2 pR \ 
2dR/dr ) 



dM 
dr 



M 



KC 2 pR 2 

AdR/dr 



dR 
dr 



(18) 



Matter density can be expressed in terms of n(z) - the observed number 
density of sources in the redshift space per steradian per unit redshift interval. 
Thus, the number of sources observed in a given redshift interval and solid angle 
is n dfl dz and the total rest mass between z and z + dz is 



M. = Anfhndz, 



(19) 



where m(z) is the average mass per source. On the other hand the total rest 
mass between r and r + dr is 



M = pd 3 V 



^4nR 2 R, r , 
P ; =dr, 
+ 



(20) 



where d 3 U is the proper volume on a constant time slice, evaluated on the null 
cone. Hence by (jT9" ]h ([20 ]) and ([12"]) 

.dz 



R p = ran 



dr 



(21) 



Finally, to find r(z), the r.h.s. of pH|) must be expressed in terms of R and n(z). 
Differentiating @ with respect to r and substituting the result in the r.h.s. of 
(fTT)]) , we obtain 



R,tr 



VT+2E R, t 



M, r 



M 



RVTT2E R 2 



+ (VT+2E),r 



(22) 



The derivative of \>T+~2E follows from (fTTJ). Then, replacing M, r by (O and 
using (|16L the above equation can be written as 



1 R, 



Now, from (jTUJ) : 



di?dz d 2 R 



Applying 



dR 
dr 



dr dr dr 2 

_ dRdz 
dz dr ' 



- kc pR 



(1 + *) 



d 2 fl 
dr 2 



di? 
dr 



--Kc 2 pR(l + z) 



(23) 



(24) 



d 2 R dR d 2 



dr 2 dz dr 2 
and integrating with respect to r, yields 



d 2 fl /dz N 2 
dz 2 \ dr 



_d 

dr 



dz dii ^ _^ 
dr dz 



1 



dr 



dr 



- kc^ p( z) R (z) ( 1 + z) — dz 
''' dz 



(25) 



Using the origin conditions [(dz/dr)(di?/dz)]o = [di?/dr)]o = 1, z(0) = and 
(|2T|) the above can be rearranged to obtain 



dz 
dr 



di? 

dz 



(1 + *) 



X < 1 KC 



m(z)n(z) 



(1 + z)dz 



(26) 
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Finally: 



r{z) = 



dF< 1 + *> 



rh(z)n(z) 



(1 + z) dz 



dz 



(27) 



3.1.2 The algorithm 

In order to specify the model, we proceed in the following way: 

1. The model is defined by two functions on the past null cone: the an- 
gular diameter distance, Da(z), and the mass density in rcdshift space, 
m(z)n(z). We assume that these functions are the same as in the ACDM 
model: 

Da{z) = — f dZ ' (28) 

37T 2 dr 

m {z)n{z) = nJ-^{l + zfD 2 A ^-, (29) 
8ttG dz 

where fl m = (8nG)/{3H$)p . 

2. Using the MHE algorithm we find r(z) by solving (f2~T|) . 

3. We numerically invert this relation to find z(r) and solve (fl~8f to find 
M(r). 

4. The function E(r) is found by solving (fT7|) . 

5. Once E and M are known, we find rj and then i^, by solving the appro- 
priate relations ([TJ) . 

6. Since in (PU) the term (1 - 2M/R)/(dR/dz) becomes 0/0 at the apparent 
horizon, the computer produces inaccurate results in the vicinity. To 
overcome this we apply the procedure described in Sec. 13.1.41 



3.1.3 The results 

The algorithm described in the previous section allows us to find an L-T model 
from a given (Da(z), m(z)n(z)) set of data. 

The free functions of the L - T model, E, ts, and M are shown in Figs. [5J 
[21 and 2] respectively. 

As can be seen, there is a problem with numerical integration for E and 
ts around r = 2.9 Gpc. The problem is related to (|17p where the term 
(1 — 2M / R) / (dR/dz) becomes 0/0 at the apparent horizon. Because of this, 
the computer produces inaccurate results in the vicinity. One solution to this 
problem was proposed by Lu & Hellaby (2007) who performed series expansions 
of R(z), n(z), dr/dz, M(z) and E(z) around the apparent horizon. However, 
this method leads either to jumps in one of these functions, say E(z), or to 
lower accuracy of the algorithm (Lu & Hellaby 2007). Therefore, we propose a 
different, much simpler approach. Namely, we fit polynomials to E(r) and M(r) 
and then we recalculate the area distance and the redshift-space mass-density 
as functions of rcdshift to check the accuracy of our approximations. As we will 
see, this method leads to results that from the observational point of view (every 
observation is accompanied with an error) are indistinguishable from those of 
the ACDM model. 



9 




0.5 1 1.5 2 2.5 3 3.5 4 
r [Gpc] 



Figure 2: The function E(r) of the L-T model defined by the {Da, m(z)n(z)) 
set of data (see Sec. I3.1.1I for details). There is a problem with numerical 
integration around r = 2.9 Gpc. The problem comes from (|17[) where the term 
(1 — 2M/R)/(dR/dz) becomes 0/0 at the apparent horizon. 




0.5 1 1.5 2 2.5 3 3.5 4 
r [Gpc] 



Figure 3: The function ts(r) of the L-T model defined by the (D^, m(z)ro(z)) 
set of data (see Sec. I3.1.1l for details). There is a cusp around r = 2.9 Gpc. The 
cusp follows from an unstable behaviour of E(r) around r — 2.9 Gpc. 
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Figure 4: The function M(r) of the L-T model defined by the [Da, m(z)n{z)) 
set of data (see Sec. I3.1.1I for details). In this case there is no problem around 
r = 2.9 Gpc. 



3.1.4 Dealing with the apparent horizon 

To overcome the numerical problem of indeterminacy of E(r) at the apparent 
horizon, we fit polynomials to the obtained M(r) and E(r). The most obvious 
choice would be polynomials in the variable r. However, in numerical experi- 
ments we noticed that much better results are obtained when we approximate M 
and E by polynomials in R(to,r), where to is the present instant. The explicit 
forms of the fitted functions are 

71=8 

M(R(t ,r)) = J2 M n^, (30) 

n— 3 

where £ = R(t , r)/lGpc, (M 3 , M 4 , M 5 , M 6 , M 7 , M 8 ) = (8.142244 x 1(T 3 kpc, 
kpc, 1.32458 x 1(T 4 kpc, -3.79 x 1(T 5 kpc, 3.23834 x 1CT 6 kpc, -9.8233 x 1CT 8 , 
kpc), and 

n— 6 

E(R(t Q ,r)) = J2 E nt n , (31) 

n=2 

where (E 2 , E 3 , E A , E 5 , E 6 ) = (1.9475 x 1CT 2 , 4.28698 x 10~ 3 , 6.50383 x 1(T 4 , 
-3.66095 x 10" 5 , 8.78679 x 10~ 7 ). 

The profiles of these functions together with those numerically derived are 
presented in Figs. [S] and |B1 We then use these functions as initial conditions 
and solve the null geodesic equations. First, we invert (|10p and (fTTj) to get 
the equations for dt/dz and dr/dz, to derive the pair (t,r) for a given redshift; 
simultaneously we solve (0) to get R, R, r , and R,t r - 

The results are shown in Figs. [7HHI The angular diameter distance is re- 
covered very accurately, while the redshift-space mass-density less so, but still 
up to z = 4 it does not differ by more than 6% from the redshift-space mass- 
density in the ACDM model — which is far less than the expected observational 
uncertainty. In addition we calculate the prediction for H(z), and we compare 
it to the estimations of the expansion rate by Simon et al. (2005). Since these 
are based on the observed age of the oldest stars, and H(z) follows from dt/dz, 
thus, as seen from (TIT]) . H(z) = R, tT jR, r . The results are presented in Fig. 
[9] As seen, the L-T model does not deviate from the ACDM model by more 
than 5%. These differences in m(z)n(z) and H(z) are caused by two factors: a) 
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Figure 5: The function £ as a function of the current areal radius for the 
(DA,m(z)n(z)) set of data - the comparison of results obtained for E as given 
by (jTTJ) (solid line) and for the approximation to E as given by (|3ip (dashed 
line). See Sec. 13.1.41 for details. The inset presents E/R 2 , the quantity that is 
constant in the FLRW limit. 



the equations (|30p and (I31[) are just approximations, b) numerical errors in the 
vicinity of the apparent horizon bias the solution of (|17|) for r > tah (where tah 
is the position of the apparent horizon). In principle, however, it is possible to 
construct the A = L-T model that matches the ACDM model. Finally, as seen 
from Fig. 1101 the current density profile does not exhibit a giant void shape. 
Instead, it suggests that the universe smoothed out around us with respect to 
directions is overdense in our vicinity up to Gpc-scales. As a consequence of our 
numerical procedure the value of density at the centre at the present time to is 
the same as the present density in the ACDM model. 



3.2 The model defined by D A (z) and H(z) 
3.2.1 The algorithm 

The algorithm used to find the L-T model consists of the following steps: 

1. The model is defined by two functions on the past null cone: the angular 
diameter distance Da(z) and the Hubble function H(z). We assume that 
these functions are the same as in the ACDM model - Da(z) is given by 
((28)) and 

H{z) = H Qy /Sl m (\ + z)* + l-n m . (32) 

2. We choose r so that ([T2l is satisfied on the past light cone of the present- 
day observer. Then using H(z) = R,t r /R, r , (fTU)) becomes 

dr 1 c 
dz = l + zH{z)' 

which for (|3"2"j) can be integrated to 
2c 1 



(34) 



(33) 
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Figure 6: The function M as a function of the current areal radius for the 
(Da, m(z)n(z)) set of data - the comparison of results obtained for M as given 
by (ITSl) (solid line) and for the approximation to M as given by (|3TJ|) (dashed 
line). See Sec. I3.1.41 for details. The inset presents M/R 3 , the quantity that is 
constant in the FLRW limit. 
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Figure 7: The angular diameter distance as a function of redshift; comparison of 
the results for the ACDM and L-T models. The inset presents the estimations 
of Da based on the type-la supernova measurements taken from the Union data 
set (Kowalski et al. 2008). See Sec. ETPl 
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Figure 8: The redshift-space mass-density as a function of redshift. The differ- 
ence between the LT and ACDM models is less than 6%. See Sec. 13.1.41 
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Figure 10: The ratio p/po of density to density at the origin. Solid lines present 
density profiles at constant t (from top to bottom: t = now, t = now — 2 x 10 9 y, 
t = now — 5 x 10 9 y). Dashed lines present density profiles at the hypersurface of 
constant age of the Universe (from top to bottom: r = now, r = now — 2 x 10 9 
y, r = now — 5 x 10 9 y). For each curve, the po is taken at the value of t or, 
respectively, r that identifies the curve. Note: all the graphs, for both foliations, 
use a comoving radial coordinate, r := R(to,r), which means that points with 
the same horizontal coordinate in the graph correspond to the same matter 
particle at all times. Only for the uppermost solid graph is R(ta, r) equal to the 
actual area-distance from the centre. 

3. Using we find p(z) and solve (|2"Tj) for n(z). 

4. We solve <|T5J) to find M(r). 

5. The function E{r) is found by solving (|17p. 

6. Once E and M are known we find r\ and then ts by solving the appropriate 
relations ©-©. 

7. As before, because of the 0/0 term in (TTT)) at the apparent horizon, we 
employ the procedure described in Sec. 13.2.31 

3.2.2 The results 

The results for E, t# and M are given in Figs. [TT] - [T3] As in Sec. 13.1. 3[ the 
functions E(r) and tg(r) evaluated by this algorithm behave unnaturally close 
to the apparent horizon, see Figs. [TT] and [T^] As before, this is caused by the 
fact that in flTTl) one has to deal with 0/0. We overcome this problem by once 
again fitting polynomials to these functions. 

3.2.3 Dealing with the apparent horizon 

The explicit forms of the fitted functions are: 




(35) 



n—3 
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Figure 11: The function E(r) of the L-T model defined by the (Da, H(zj) set 
of data (see Sec. I3.2.1l for details). Around r = 2.9 Gpc there is a problem with 
the numerical algorithm. 
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Figure 12: The function ts(r) of the L-T model defined by the (Da, H(z)) set 
of data (see Sec. 13.2.11 for details). 
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Figure 13: The function M(r) of the L-T model defined by the (Da, H(z)) set 
of data (see Sec. 13.2.11 for details). In this case there is no problem around 
r = 2.9 Gpc. 



where (M 3 , M 4 , M 5 , M 6 , M 7 ) = (8.2 x 10 3 kpc, -1.1948 kpc, 1.07521 x 10 2 kpc, 
-24.1385 kpc, 1.15743 kpc), and 

n— 6 

E(R(t 0) r)) = J2 E nt n , (36) 

n=2 

where {E 2 ,E 3 , E 4 , E 5 .E e ) = (1.7324x 1(T 2 , -2.5725x 10~ 3 , 1. 14925 x 10~ 4 , 2.87776x 
10- 5 , -1.90389 x 10- 6 ). 

The profiles of these functions together with the numerically derived ones 
are shown in Figs. H4l and [T5l We then use these functions as initial conditions 
and solve the null geodesic equations. The results are presented in Figs. [TBI and 
1171 Both the angular diameter distance and the expansion rate as functions of 
the redshift are recovered very accurately. From the observational perspective 
these two models are indistinguishable. In addition we present the m(z)n(z) 
plot. As seen from Fig. [TO] it also gives quite an accurate fit, with a deviation 
from the ACDM model being less than 2.5%. Finally, as seen from Fig. [TO] the 
current density profile has a similar shape as in Sec. 13.1.41 and this is not a 
giant void. 

4 Discussion and conclusion 

Contrary to what is commonly claimed, L-T models with a giant void do not 
reproduce the main features of the ACDM model. These types of models just fit 
cosmological observations, with a priori constraints imposed on the L-T models. 
Indeed, we have found the L-T models that mimic some of the observational 
features of the ACDM model and exhibit no giant void, but rather a giant hump. 

It is clear from the energy density profiles shown in Figs. [TO] and [TO] that 
the L-T models we have reconstructed have a large overdensity, when viewed 
over large enough scales. However, as we have made clear throughout, these 
profiles are the result of reproducing observables that match ACDM predictions, 
and not from fitting to any real data. Luminosity distances from supernovae 
observations, and number counts from galaxy surveys do not currently extend 
much beyond z ~ 1.5 and z ~ 0.5, respectively. As such, using real observables, 
it would only currently be possible to perform a reconstruction of a limited part 
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Figure 14: The function E as a function of the current areal radius for the 
(Da,H(z)) set of data - the comparison of results obtained for E as given by 
(ITT)) (solid line) and for the approximation to E as given by (|5fJ)l (dashed line) . 
See Sec. l3.2.3l for details. The inset presents E/R 2 , the quantity that is constant 
in the FLRW limit. 
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Figure 15: The function M as a function of the current areal radius for the 
(Da, H(z)) set of data - the comparison of results obtained for M as given 
by (|T51) (solid line) and for the approximation to M as given by (dashed 
line). See Sec. 13.2.31 for details. The inset presents M/R 3 , the quantity that is 
constant in the FLRW limit. 
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Figure 16: The angular diameter distance as a function of z. Sec Sec. 13.2.51 for 
details. 
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Figure 17: The function H(z). See Sec. 13.2. 31 for details. For comparison the 
measurements of H (z) (Simon et al. 2005) are also shown. 
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Figure 18: The redshift-space mass-density as a function of redshift for the 
LT model considered in Sec. 13.2.31 and in the ACDM model. The difference 
between the two models is less than 2.5%. Compare Fig, [TJ the consistency 
achieved here is much better than in a giant void model. 
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Figure 19: The ratio p/po of density to density at the origin. Solid lines present 
density profiles at constant t (from top to bottom: t = now, t = now — 2 x 10 9 y, 
t = now — 5 x 10 9 y). Dashed lines present density profiles at the hypersurface of 
constant age of the Universe (from top to bottom: r = now, t — now — 2 x 10 9 
y, t = now — 5 x 10 9 y). For each curve, the po is taken at the value of t 
or, respectively, r that identifies the curve. The remark about the horizontal 
coordinate in Fig. [TU] applies also here. 
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of the full structure we have found here, and even then only up to the degree 
allowed by the errors associated with these quantities. If one were to attempt 
such a reconstruction, it appears from Figs. 1101 and [T^l that one may, in fact, 
reconstruct a local energy density profile that is an increasing function of r, 
and not decreasing, due to the limited extent of these observations in z (this 
is particularly true in the equal age foliation). Our interpretation of a giant 
hump should therefore be understood as corresponding to an extrapolation of 
observable quantities under the expectation that they will follow ACDM, rather 
than being directly implied by any currently known observations themselves. 

Recently, some astrophysicists have begun to take seriously the cosmological 
implications of the L-T model. This model, although still simplejf] is quite pow- 
erful and exhibits some features of general relativistic dynamics, like arbitrary 
functions in the initial data. 

As we said earlier in this paper, the belief that an L-T model fitted to 
supernova la observations necessarily implies the existence of a giant void with 
us at the centre was created by arbitrarily limiting the generality of the model. 
With its free functions fitted to ACDM features rather than to expectations, the 
giant void does not necessarily follow. Rather, one alternative is that the graph 
of the density smoothed out over angles around us has the shape of a shallow 
and wide valley on top of a giant hump. 

This giant hump may be a feature of the particular L-T model that we ended 
up with. Variations in the ACDM parameters would modify the details, and the 
observational constraints we considered are not very tight. Future calculations 
with other constraints may favor a still different profile at t — now. Hence 
we do not wish our paper to become a starting point of a new paradigm in 
observational cosmology, aimed at detecting the hump. Before this happens, we 
must decide at the theoretical level whether the hump is a necessary implication 
of L-T models properly fitted to other observations. 

It must be stressed that this hump is not directly observable. It exists in the 
space t = now, of events simultaneous with our present instant in the cosmolog- 
ical synchronisation, i.e. it is in a space-like relation to us. This is also the case 
of the giant void (see, e.g., Fig. 4 of Alnes et al. 2006; Fig. 1 of Garci'a-Bellido 
and Haugb0lle 2008a; Figs. 4 and 6 of Yoo et al. 2008). However, an observa- 
tional test of the giant void is easier to complete. The reason is that the models 
considered in this paper have a redshift-space mass-density almost the same as 
in the ACDM model (see Figs. 151 and fT5|) and as can be seen from Fig. [201 their 
p(z) scales with the redshift in almost the same manner as in the FLRW model, 
i.e. ~ (1 + z) 3 . This is not the case of the giant voicjf] which does not reproduce 
these features on the past light cone (see Figs. IT1 and USD")) . Thus, unlike the giant 
void, the giant hump is not observable in p(z) or in the number count data. 

What is the cause of this difference between the density distribution on 
our past light cone and in the t = now space? It is the oft-forgotten basic 
feature of the L-T model (and in fact of all inhomogeneous models, also those 
not yet known explicitly as solutions of Einstein's equations): on any initial 
data hypersurface, whether it is a light cone or a t = constant space, the density 

4 From the computational point of view, the difference between the Friedmann and L-T 
models is quite trivial. The equation ||2} that governs the evolution of the L-T model is 
exactly the same as in the Friedmann model; it is still an ordinary differential equation in the 
time- variable t. The only difference is that the function R(t) obeying J2]l depends on one more 
variable, the radial coordinate r. This r enters only as a parameter, and then automatically 
all the integration 'constants' that appear while solving J2j are no longer constant, but are 
functions of r. Sophistication comes at the level of interpreting the solutions - however, this 
is no longer mathematics, but astrophysics. 

5 The giant void used here is Bolejko & Wyithe (2009)'s model with radius of 2.96 Gpc and 
density contrast of 4.05. The redshift-space mass-density for this model is presented in Fig. 
1. 
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Figure 20: The density distribution as a function of redshift for the giant void 
(GV) L-T model, for the models of Sec. 15X1 and l5T5l (GH) and for the FLRW 
model, where p(z) = po(l + z) 3 (dashed line). The parameters of the GV model 
are: radius 2.96 Gpc and density contrast 4.05. The same model is presented 
in Fig. 1; for details of how the model is specified and null gedesics solved see 
Bolejko & Wyithe 2009. 



and velocity distributions are two algebraically independent functions of position. 
Thus the density on a later hypersurface may be quite different, since it depends 
on both initial functions. Whatever initial density distribution we observe can be 
completely transformed by the velocity distribution. For example, as predicted 
by Mustapha and Hellaby (2001) and explicitly demonstrated by Krasihski and 
Hellaby (2004), any initial condensation can evolve into a void and vice versa. 
In FLRW models, there are no physical functions of position, and all worldlines 
evolve together. Thus, while dealing with an L-T (or any inhomogeneous) 
model, one must be cautious when applying the Robertson- Walker-inspired 
prejudices and expectations. 

The use of oversimplified L-T models can create another false idea and false 
expectation. The false idea is that there is an opposition between the ACDM 
model, belonging to the FLRW class, and the L-T model or in general, inho- 
mogeneous models: either one or the other could be 'correct', but not both. 
This putative opposition can then give rise to the expectation that more, and 
more detailed, observations will be able to tell us which one to reject. In truth, 
there is no opposition. The inhomogeneous models, like for example the L-T 
model with its two arbitrary functions of one variable are huge, compared to 
FLRW, families of models that include the Friedmann models as a very simple 
subcase. The fact, demonstrated in several papers already (see Celerier 2007, 
for a review), that even a A — L-T model can imitate A ^ in an FLRW 
model, additionally attests to the flexibility and power of the L-T model. Thus, 
if the Friedmann models, ACDM among them, are considered good enough for 
cosmology, then the L-T models can only be better: they constitute an exact 
perturbation of the Friedmann background, and can reproduce the latter as a 
limit with an arbitrary precision. While future observations, for example the 
kSZ effect (Garcia-Bellido & Haugb0lle 2008b) or the growth of linear struc- 
ture (Clarkson et al. 2009)) will provide a sufficient insight to test particular 
configurations (like for example a giant void model), we will never be able to 
reject inhomogeneous models. After all, the Universe as it is, is inhomogeneous. 
Nowadays we use homogeneous models just for simplicity, and although they 
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have worked well so far, in future they will certainly be replaced by more so- 
phisticated models, either by exact solutions, or what is more probable in light 
of increasing computation power of computers, by numerical simulations. 

When considering models that go beyond the FLRW approximation, one 
may ask either 'what limitations on the arbitrary functions in the models do 
our observations impose', or 'which model best describes a given situation: a 
homogeneous FLRW model or an inhomogeneous one?' The latter of these 
questions has often been asked in the context of comparing L-T models without 
A to FLRW models with A, and is of much interest for understanding the ne- 
cessity of introducing A into the observer's standard cosmological model. Such 
hypothesis testing questions are often posed in cosmology, but are difficult to 
address in the current context as they require artificially limiting the generality 
of the models in question (in order to have a finite number of parameters, so 
that the test can be performed) . Given the lack of motivation for exactly how 
to perform such a limitation, one is then left in the undesirable circumstance of 
(potentially) dismissing particular L-T models, while being left with an infinite 
number of remaining L-T models to evaluate. We therefore prefer to consider 
the former question. In order to reasonably answer this for the L-T model, a 
general framework for interpreting observations in the L-T geometrical back- 
ground should be created (and in the future it should be transformed into a 
framework for interpreting the observations in a still more general, or the most 
general geometrical background). Such a program is still in its infancy, but is 
being actually developed by C. Hellaby and coworkers under the name 'Metric 
of the Cosmos' (Lu and Hellaby, 2007, McClure and Hellaby, 2008, Hellaby and 
Alfedeel, 2009). 
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